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Abstract 

We have investigated the magnetic properties of a piezoelectric actuator /ferromagnetic semicon- 
ductor hybrid structure. Using a GaMnAs epilayer as the ferromagnetic semiconductor and ap- 
plying the piezo-stress along its [110] direction, we quantify the magnetic anisotropy as a function 
of the voltage V p applied to the piezoelectric actuator using anisotropic magnetoresistance tech- 
niques. We find that the easy axis of the strain- induced uniaxial magnetic anisotropy contribution 
can be inverted from the [110] to the [HO] direction via the application of appropriate voltages 
V p . At T = 5 K the magnetoelastic term is a minor contribution to the magnetic anisotropy. 
Nevertheless, we show that the switching fields of p{^qH) loops are shifted as a function of V p at 
this temperature. At 50 K — where the magnetoelastic term dominates the magnetic anisotropy 
— we are able to tune the magnetization orientation by about 70° solely by means of the electri- 
cal voltage V p applied. Furthermore, we derive the magnetostrictive constant Am as a function 
of temperature and find values consistent with earlier results. We argue that the piezo-voltage 
control of magnetization orientation is directly transferable to other ferromagnetic/piezoelectric 
hybrid structures, paving the way to innovative multifunctional device concepts. As an example, 
we demonstrate piezo-voltage induced irreversible magnetization switching at T = 40 K, which 
constitutes the basic principle of a nonvolatile memory element. 

PACS numbers: 75.30.Gw, 75.50.Pp, 75.47.Pq, 77.65.-j 
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I. INTRODUCTION 



Ferromagnetic semiconductors unite the long-range magnetic ordering characteristic of 
ferromagnets with the versatile properties of conventional semiconductors. This class of mul- 
tifunctional materials therefore is very attractive from a fundamental physics point of view. 
Moreover, novel spin-electronic devices can be realized using ferromagnetic semiconductors, 
in which the electronic functionality is directly linked to magnetic properties such as the 
magnetization orientation For device applications it is attractive to 

control the magnetic properties via a non-magnetic parameter. In this paper, we investigate 
the approach to strain a dilute magnetic semiconductor (DMS) lattice by means of an ex- 
ternal piezoelectric actuator in order to manipulate its magnetic properties. Using GaMnAs 
as the DMS, we discuss in detail the piezo-voltage control of the magnetic anisotropy, which 
can be utilized to reversibly manipulate magnetization orientation as well as to irreversibly 
switch the magnetization between two magnetic easy axis orientations. The experiments 
discussed here pave the way for the three dimensional control of spins via a non-magnetic 
control parameter. 

We begin with a short introduction into GaMnAs, before reviewing different approaches 
to influence the magnetic properties of magnetic semiconductors via non-magnetic param- 
eters reported in literature. Gai-^Mn^As is the prototype ferromagnetic semiconductor. 
Although its Curie temperature Tq < 170 K is still below room temperature, it has been 
investigated vigorously in the last decade The effect of the Mn incorporation is twofold: 
On the one hand, the half-filled <i-shells of the Mn 2+ ions provide localized magnetic moments 
with spin S = 5/2. On the other hand, Mn also is a shallow acceptor in GaAs, so that ferro- 
magnetic GaMnAs doped with a few percent of Mn is degenerately p-type. These itinerant 
holes mediate the ferromagnetic exchange between the localized Mn moments , linking; 
the magnetic properties of GaMnAs to the GaAs valence band structure sj. The strong p-d 
coupling between the delocalized holes and the localized electrons in the Mn ci-shell results in 
large magnetoresistive effects in GaMnAs. For example, the anisotropic magnetoresistance 
(AMR) has a "giant" transverse (planar Hall) component [u|. The anomalous Hall effect 
due to the magnetization in GaMnAs is routinely exploited for magnetic characterization 



purposes 



11] . Moreover, the spin polarization of the itinerant holes is large, with values 



of up to 85% reported in literature [l2, [ijj]. Since the magnetic domains extend over 
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100 microns or more [151 ] . GaMnAs is ideally suited as a spin-injecting or spin-polarizing 
contact in spintronic devices. Indeed, novel functionalities or physical effects have already 
been demonstrated in spin-electronic devices based on GaMnAs, such as the emission of 
circularly polarized light from a so-called spin light emitting diode {12I . [if] ], the controlled 



motion of domain walls via current pulses 17|, ll8( , the tunneling anisotropic magnetoresis 



tance in GaMnAs/insulator /normal metal structures [19j], and non- volatile memory device 
concepts Q, 7 , Q] . 

Gai^Mn^As layers grown on GaAs are compressively strained, leading to a first-order 
uniaxial magnetic anisotropy so that the film plane is an easy magnetic plane for the typical 



22] . A smaller cubic (biaxial) anisotropy 



range of Mn concentrations (0.01 < x < 0.1) 
yields two approximately orthogonal easy axes within the film plane. Furthermore, a small 
in-plane uniaxial anisotropy along [110] is present in most GaMnAs films on GaAs 23 
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29. 
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31] , the microscopic origin of which is still controversially discussed. In 



GaMnAs films with tensile strain, in contrast, the film plane is magnetically hard 20 



so 



that the magnetization aligns along the film normal at low magnetic fields. Furthermore, the 
magnetic anisotropy in GaMnAs qualitatively changes with the temperature. In comparison 
to other ferromagnetic metals, demagnetization effects only play a minor role in GaMnAs, 
as the saturation magnetization is small due to the small concentration of magnetic ions. 

One of the most intriguing properties of DMS is the strong dependence of their magnetic 
pro per ties on non-magnetic parameters, such as electric field 321. 133 ], light irradiation 34, 



35 



36] . temperature [24j, 



26 



37] . dopant density [38, 



39] . strain 



in 3 



21 



22 



40], or pressure 



4=11 ] . This allows for various control schemes of the magnetic properties of DMS, which we 



shortly discuss in the following. 



Ohno et al. 



321 ] demonstrated that an electric field-control of ferromagnetism is possible 



in InMnAs samples covered with a metallic gate. Upon the application of an appropriate 
voltage to the gate electrode, the density of itinerant holes is altered, which allows to switch 
on or off the ferromagnetic exchange. Electrically assisted magnetization reversal, as well 
as electrical demagnetization have already been demonstrated in devices utilizing this effect 



331 ] . While this electric field-control of ferromagnetism appears very attractive for spintronic 



devices, the high carrier density in degeneratly doped DMS and the resulting screening effects 
limit the thickness of magnetic semiconductor films controllable by electric fields to a few 
nm at most, imposing severe constraints on the device geometries. 
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Another approach to vary the hole concentration and thus the magnetic properties relies 
on the introduction of additional non-magnetic dopant species into GaMnAs. By co-doping 
GaMnAs with Be, the relative strength of the cubic and uniaxial magnetic anisotropy contri- 



butions can be altered 



The incorporation of up to 17% of Al into GaMnAs significa ntly 

anc 



reduces the density of free holes, eventually leading to a magnetically hard film p 



The incorporation of hydrogen also allows to contro 



the hole density in GaMnAs 



39]. 



421. The 



hydrogen atom forms a complex with the Mn ion [431 . |44j and thereby passivates its acceptor 
function, while the localized 3c? electrons stay unaffected. The hydrogen incorporation is sta- 
ble at ambient temperatures, but can be reversed by a mild thermal as well as a pulsed-laser 



annealing 45j, |46| . However, while the exchange coupling in GaMnAs films with thicknesses 



of several hundred nm can be switched on or off via co-doping or hydrogenation, these tech- 
niques do not allow to manipulate the magnetization orientation while operating a spintronic 
device. 

For the realization of DMS tunable during device operation, the dependence of the mag- 
netic properties on the crystal lattice appears much more promising. Csontos et al. showed 
that the ferromagnetic exchange interaction in InMnSb can be controlled via the application 

n 

of hydrostatic pressure [41]. As discussed above, the magnetic anisotropy in GaMnAs was 



found to qualitatively change function of strain 



20 



21 



22 



4=01 ] . Recently also the con- 



trol of the mangnetic anisotropy via anisotropic strain relaxation in patterened structures 
has been reported j^]. These experimental observations are explained by the established 
theoretical model for the ferromagnetic exchange 9]: The strain alters the symmetry of the 
hole wave function in the GaAs valence band, which due to the hole-mediation of ferro- 
magnetic exchange has a strong influence on the magnetic anisotropy. This picture also is 
compatible with the recent measurements of the magneto-elastic constants of GaMnAs in 
micro-mechanical devices 37 1. 



In the following we will investigate the control of the magnetic anisotropy of GaMnAs 
via the application of an external stress by a piezoelectric actuator. This approach has 
already been shown to allow in-situ manipulation of magnetic anisotropy in piezoelectric ac- 



tuator /ferromagnet hybrid samples at room temperature 48|, |49j. We have recently demon- 



strated that this approach can be transfered to piezoelectric actuator/GaMnAs hybrids 50], 



as also reported by 



51 



521 ] . Here we compile our extensive experimental findings and show 



that they can be understood consistently using a free energy model. In particular, we make 
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use of the strong temperature dependence of the magnetic anisotropy in GaMnAs to tune 
the free energy surface so that the effect of strain on the magnetic properties is maximal. 

This paper is organized as follows: After a review of the magnetic free energy approach 
in Sec. [IT] we summarize the experimental procedures employed and the properties of the 
piezoelectric actuator used in Sec. II III and IIV[ respectively. The method to determine the 
magnetic anisotropy from magnetotransport experiments is described in Sec. [V] At low 
temperatures (T m 5 K) the piezo- voltage induced changes are small compared to the domi- 
nating cubic contribution to magnetic anisotropy in GaMnAs. In this temperature range we 
demonstrate the piezo-voltage control of the magnetoresistance switching fields in external 
magnetic field sweeps (Sec. I VI p . In contrast, at T « 50 K the magnetoelastic contribution 
dominates magnetic anisotropy (Sec. IVIII) . which allows to continuously and reversibly ro- 
tate magnetization orientation (Sec. IVIIII) . From a comparison of the piezo-induced strain 
(Sec. lIVp and the corresponding magnetoelastic contribution to magnetic anisotropy we are 
able to derive the temperature dependence of the magnetostrictive constant (Sec. HXl) . Fi- 
nally, we demonstrate a piezo-voltage induced irreversible magnetization switching (Sec. IXl) 
and discuss possible applications of the investigated effects (Sec. IXB . 



II. FREE ENERGY APPROACH 



To describe the magnetic anisotropy and thus magnetization orientation in our sample 



we adopt a free energy approach. To seconc 
to the saturation magnetization is given by 



order, the free energy in GaMnAs normalized 



21 



53] 



F/M = -fi H(h ■ m) + B c \\(mi + mj) + B 110 (V p )(j ■ m) 2 + B 010 m 2 y = 
= -/i H cos(/3 - a) + B cll [cos 4 (/3 + 45°) + cos 4 (/3 - 45°)] + 

+B U0 (V P ) cos 2 + £? io cos 2 (/3 - 45°), (1) 

where we have assumed that the magnetization is aligned within the film plane, /io = 
4tt x 1(T 7 is the vacuum induction, m,- are the direction cosines of the magnetization 

Am ' 1 ° 

relative to the cubic axes, B c \\ is the in-plane cubic anisotropy parameter, and Bno(V p ) and 
£>oio are the uniaxial anisotropy parameters. The orientation of the external magnetic field 
H and the magnetization M are given by the unit vectors h and m, with the corresponding 
angles a and (3 with respect to the current direction j as shown in Fig. [3j The individual 
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FIG. 1: (color online) (a) Schematic illustration of a 90° magnetization switching for a GaMnAs 
film affixed onto a piezoelectric actuator with the magnetic easy axis [100] parallel to the main 
expansion direction of the actuator, (b) Corresponding free energy surface. The free energy for 
the magnetization oriented along a specific direction is given by the distance from the center of the 
coordinate system to the surface. In this example the magnetization is oriented along [100] . (c) A 
positive voltage applied to the actuator leads to an elongation of the actuator/GaMnAs layer, (d) 
This strain induces a uniaxial anisotropy with a magnetic hard axis in elongation direction, which 
leads to a change in the relative strength of the in-plane magnetic easy axes [Fig. [2{a)] and thus a 
90° magnetization switching into the global minimum at [010] . 

terms in ([1]) describe the Zeeman interaction, the in-plane cubic and two uniaxial magnetic 
anisotropy contributions along [110] and [010] crystallographic directions. The third term 
[i?iio(Vp)(j ■ m) 2 ], the magnetoelastic contribution to the magnetic anisotropy, accounts for 
the changes in magnetic anisotropy due to a voltage V p applied to the piezoelectric actuator 
in the particular stress geometry discussed below. The fact that the piezo-voltage induced 
effect can be described by a uniaxial anisotropy is derived in the Appendix. Note that the 



i 2K U 

anisotropy fields used in Ref. |54|, [55|] are given by = —4B C \\ for the in-plane cubic field 
and — = 2Bi for the uniaxial fields. 

Figure [TJb) illustrates a typical free energy surface of a compressively strained GaMnAs 
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FIG. 2: (color online) Cuts through the free energy surface within the film plane for two different 
strain configurations. The black dashed curves illustrate the in-plane cubic anisotropy contribution 
with magnetic easy axes along (100). (a) For the main elongation direction of the piezoelectric 
actuator aligned along a magnetic easy axis the strain-induced uniaxial anisotropy leads to a change 
of the relative strengths of the magnetic easy axes (red full curve), (b) In contrast, for the main 
elongation direction of the piezoelectric actuator aligned along a locally hard axis of the in-plane 
cubic anisotropy, e.g. [110], the strain-induced uniaxial anisotropy leads to a rotation and a change 
of the relative orientations of the magnetic easy axes (red full curve). 

film. The free energy for the magnetization oriented along a specific direction is given 
by the distance from the center of the coordinate system to the surface. The magnetic 
anisotropy in Fig. [T](b) is dominated by a uniaxial anisotropy perpendicular to the film 
plane (the film plane is indicated by the orange plane). Consequently, the magnetization 
will be oriented within the film plane at low external fields in good approximation. In the 
film plane magnetic anisotropy mainly is determined by the cubic contribution with easy 
axes along (100) at liquid He temperatures. 

To achieve a piezo-control of the magnetization orientation, the GaMnAs film is cemented 
onto the piezoelectric actuator (Sec. IIIII) . Regarding the relative orientation of the GaMnAs 
film and the main elongation direction of the piezoelectric actuator there are two qualita- 
tively different approaches: (i) To achieve an irreversible switching of the magnetization 
from one minimum of cubic anisotropy into another, the GaMnAs film has to be affixed 
onto the piezoelectric actuator with an easy axis parallel to the main expansion direction of 
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FIG. 3: (color online) Illustration of the relative alignment of the Hall bar, the GaMnAs thin film, 
and the main expansion direction of the piezoelectric actuator, h = H/H, m = M/M, j, and t 
denote the unit vectors along the orientation of the magnetic field H, the magnetization M, the 
current density j, and the direction transverse to the current, respectively. 



the actuator [Fig. QJa)] 52j|. In this way, starting from a GaMnAs film magnetized along 
the main expansion direction, the application of a positive voltage will lead to tensile strain 
and thus induce a uniaxial magnetoelastic contribution with a magnetic hard axis along the 
expansion direction [Fig. DJc)]. Thereby, the magnetization can be switched by 90° into the 
direction of the new global minimum in free energy surface [Fig. [2(d), Fig. 0(a)] . We note 
that two equivalent global minima are present at [010] and [010] . In order to lift this degen- 
eracy, e.g. a small external magnetic field should be applied in experiment, (ii) In contrast, 
to be able to continuously and reversibly rotate magnetization within the film plane, we 
here will affix the GaMnAs film onto the piezoactor with the orientation of the [110] sample 
edge parallel to the main expansion direction (Fig. [3]). The difference between these two 
approaches can be visualized via cuts through the free energy surface within the film plane 
(Fig. [2]). While in the alignment (i) [Fig. 12(a)] the relative strength of the magnetic easy 
axes is changed, the alignment (ii) [Fig. [2(b)] allows a rotation and thus a change of the 
relative orientation of the magnetic easy axes. 
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III. EXPERIMENTAL 



The 30 nm thick Gai-^Mn^As film with a Mn concentration x = 0.045 investigated 
here was grown on a (OOl)-oriented GaAs substrate by low-temperature molecular-beam 
epitaxy (LT-MBE) j^j. The sample exhibits a Curie temperature Xc ~ 85 K determined 
from superconducting quantum interference device (SQUID) magnetometry. Via optical 
lithography and wet chemical etching we patterned the film into Hall bars with the current 
direction j along the [110] crystal axis. After the etching process the GaAs substrate was 
mechanically polished to a thickness of about 100 /xm. The sample then was cemented onto 
a lead zirconate titanate (PZT) piezoelectric actuator "PSt 150/2 x3/5a" (Piezomechanik 
Miinchen) using the two-component epoxy "M-Bond 600" (Vishay Inc.) annealed for 2 h 
at 120 °C in air. As shown in Fig. [3] the main piezoelectric actuator expansion direction is 
aligned parallel to the GaAs [110] direction. 

We extracted the magnetic anisotropy from magnetotransport experiments carried out in 
a superconducting magnet cryostat. The sample was mounted on a rotatable sample stage 
that allows rotation around one axis, in such a way that the applied magnetic field always lay 
within the film plane, and that the angle a enclosed by the magnetic field orientation h and 
the current density j could be adjusted within —140° < a < 140°, as illustrated in Fig. [31 
Using a dc current density of j = 4.4 x 10 3 Acm~ 2 the resistivities pi ong along j and ptrans per- 
pendicular to j were recorded via four-point measurements. We studied the magnetic prop- 
erties of the sample using conventional magnetoresistance traces {p\ ong (fioH) , ptransi/J-oH)} 
at a fixed magnetic field orientation a, as well as via angle- dependent magnetoresistance 
(ADMR) measurements {p long (a), Ptrans(«)} at a fixed external magnetic field strength. 

IV. MEASUREMENT OF PIEZO-INDUCED STRAIN 

Application of a positive piezo-voltage results in a lattice elongation 6jj along j and a 
lattice contraction e a along t perpendicular to j (cf. Fig. [3]). To determine the temperature 
dependence of the piezo-induced lattice distortions en(Vp), i G {j, t} we used strain gauges. 
Figure @Ja) depicts the measured longitudinal (e^) and transverse (e 4f ) lattice distortions 
for piezo-voltage sweeps at 5 K and 50 K. For increasing temperature we observe an increase 
in both the total lattice distortions Ae^ = eu(V p = +200 V) - eu(V p = -200 V), i G {j,t} 
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FIG. 4: (color online) (a) Longitudinal (e,-j, black full symbols) and transverse (e^, red open sym- 
bols) lattice distortions for piezo-voltage sweeps at 5 K (squares) and 50 K (circles) determined 
via strain gauges. The data show three full voltage cycles. To allow for comparison, we plot 
£u (V p ) — en (+200 V), i € {j,t}. (b) Temperature dependence of the piezo-induced lattice distor- 
tions Ae^i = ea (Vp = +200 V) - €u(Vp = -200 V), i G {j,t}, determined from the third voltage 
cycle. 

[Fig. S(b)] and the hysteresis exhibited by the piezoelectric actuator [Fig. IU(a)]. This piezo- 
voltage induced strain causes an additional magnetoelastic contribution F mage i(V p )/M = 
Bno(Vp) cos 2 (3 to the magnetic anisotropy as derived in the Appendix. Positive (negative) 
piezo-voltage thus results in an additional uniaxial magnetic anisotropy with a magnetic 
hard (easy) axis along j-direction. 



V. DETERMINATION OF MAGNETIC ANISOTROPY FROM MAGNETO- 
TRANSPORT 

We now describe the determination of the magnetic anisotropy from magnetotransport 
measurements in which the external magnetic field H is rotated at different field strengths 
within the film plane as introduced above. Due to the dominating uniaxial anisotropy in 
growth direction caused by epitaxial strain, the magnetization in our film in good approx- 
imation is aligned within the film plane. Therefore, the longitudinal and the transverse 
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resistivities 



Piong = Po + Pi(j • m) + p 3 (j • m) = p + Pi cos (3 + p 3 cos (3 (2) 
Ptrans = /^(t ■ m) (j ■ m) = -p 7 sin(2/3) (3) 

are determined by the angle f3 between the in-plane orientation of the magnetization m and 
the current direction j. The orientation t denotes the in-plane direction perpendicular to j 
(Fig- ED- These expressions can be derived from a series expansion of the resistivity tensor in 



powers of the magnetization components m; 



53] . Following Limmer et al. 53j the resistivity 



parameters po, pi, P3, and p^ can be deduced from magnetotransport measurements in 
which a constant large external magnetic field H is rotated in the film plane. In good 
approximation, the magnetization in this case can be assumed to be oriented along the 
magnetic field direction, i.e. [3 = a in Eq. (j2J) and where a denotes the angle between 
the magnetic field orientation h and the current direction j. Equations (j2J) and ([3]) then 
yield the black dashed curves for {pi ong (a), Ptrans («)} in Fig. [51 Note that angles are defined 
in the mathematically positive sense (counterclockwise) in Fig [3] and are plotted clockwise 
in the following. Thus +90° corresponds to the [110] direction and not [l 10] as given in 
{5(|. In contrast, for measurements at smaller external magnetic fields the orientation of 
magnetization — and also pi on g and ptrans — to an increasing extent will be influenced by 
the magnetic anisotropy. Therefore, at small magnetic fields the magnetization will tend 
to remain oriented close to a magnetic easy axis (e.a., depicted as green straight lines in 
Fig. [5]) — and pi on g and ptrans will tend to remain constant — over a broad range of external 
magnetic field orientations near to a magnetic easy axis. Accordingly, abrupt changes in 
{piong(a), Ptrans(a)} indicate a nearby magnetic hard £IX1S ; clS shown by the red full curves 
in Fig. [5j Note that a similar approach to determine the in-plane magnetic anisotropy in 



56) 



GaMnAs has recently been followed by Yamada et al. 

For the angle-dependent magnetoresistance (ADMR) measurements we first aligned the 
magnetization into a well-defined initial state in an external magnetic field of PqH = 7 T 
along a = —140°. Then, we lowered the field to the measurement field and started the 
angular scan {pi on g (a) > Ptrans («)}• Figure [6] shows { pi n g ( a), Ptrans («)} measured at T = 5 K 
and V p = V for rotations within external magnetic fields pqH of 1 T, 250 mT, 100 mT, 
and 60 mT. Following the discussion above, the angular dependence at small p H already 
indicates that the in-plane anisotropy is dominated by a cubic contribution with magnetic 
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f . 180° 180" 

(a) ^ — ' — l (b) 



Plong ^ rb - UnitS > Ptrans ^ rb - UnitS > 

FIG. 5: (color online) Black dashed curves: Illustration of (a) longitudinal resistivity p\ ng( a ) an d 
(b) transverse resistivity ptrans(«) for a large external magnetic field H rotated within the film 
plane. Since in this case the magnetization is aligned in field direction, i.e. (3 = a, the resistivities 
can be calculated from Eq. (J2j) and (]3]). The parameters have been chosen with p$ and pi > 
and with p\ and p^ < 0. Red full curves: At small magnetic fields H, the magnetization tends 
to remain oriented close to a magnetic easy axis (e.a., green straight lines) and pi on g an d ptmns 
therefore tend to remain constant over a broad range of external magnetic field orientations nearby 
a magnetic easy axis. 

easy axes along (100). This is corroborated by the simulation, where the simulated curves 
(yellow dashed) are the result of a procedure in which the resistivity parameters and the 



anisotropy parameters from Eq. <^ are fitted iteratively [53J. The values for p , p 1; p 3 , and 
Pt and the anisotropy fields thus obtained are summarized in Fig. [13] below. We note that 
in this fit all resistivity parameters except for po are kept constant for the different values of 
the applied magnetic field. The change in p with magnetic field strength accounts for the 



influence of negative magnetoresistance typically observed for GaMnAs 571 ] . 



VI. ADMR AND PIEZO-STRAIN AT 5 K 

At 5 K the influence of piezo-induced strain on the ADMR scans {pi on g(aO, Ptrans(«)} dis- 
cussed in the preceding section is only marginal. This is due to the small piezo-induced 
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FIG. 6: (color online) (a) Longitudinal resistivity pi ong (a) and (b) transverse resistivity ptrans(a:) 
for a constant external magnetic field poH = 60 mT, 100 mT, 250 mT, and 1 T rotated within 
the film plane for V p = V. The black circles represent the experimental data, the yellow dashed 
curves the result of the corresponding simulation using the parameter values given in Fig. [T3j 

lattice distortion at 5 K (Fig. UJ) and the large cubic anisotropy at this temperature. 
However, the influence of the magnetoelastic contribution can be resolved in conventional 
{p\ ong (p H), Ptmnsif^oH)} magnetotransport measurements. 

Figure [7| shows the {pi on g(poH) , p tTSljns (fioH)} downsweep curves measured at 5 K for the 
external magnetic field H aligned at a = 350°. The evolution of the downsweep curve for 
V p = V from positive to negative magnetic fields can be understood as follows (compare 
the corresponding evolution of {pion g («), Ptrans(«)} m Fig. [6] at p H = 1 T where j3 = a is 
fulfilled): (A) p^H = +100 mT-> —5 mT: rotation of magnetization from (3 = 350° — > 315°, 
the orientation of the closest magnetic easy axis. (B) p H = —5 mT: first switching process 
by A/3 ~ 90° into a direction close to the magnetic easy axis at 225°. We refer to the 
corresponding field as a switching field in the following. (C) poH = —5 mT^ —40 mT: 
slight rotation of the magnetization (corresponding to a decrease of (3) within the minimum 
of free energy close to 225°. (D) p H = —40 mT: second switching process by A/3 ~ 90° 
into a direction close to the magnetic easy axis at 135°. (E) p H = —40 mT-> —100 mT: 
rotation of magnetization from (3 ~ 135° — > 170°, the orientation of the external magnetic 
field. Note that the change in {pi ong (p H), Ptrans(^o-^)} during each step (A to E) of the 
downsweep described above is in agreement with the values expected for the corresponding 
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n H (mT) n H (mT) 

FIG. 7: (color online) (a) |0i O ng (Po-^O and (b) Ptrans(Po-£0 curves measured at 5 K for the external 
magnetic field H aligned at a = 350°. An explanation of the evolution is given in the text. 

orientation of magnetization from the angle- dependent resistivity loops {pi ong (tt), Ptrans(a)} 
at high external magnetic field in Fig. El where f3 = a is fulfilled. For example, this com- 
parison explains the fact that the jumps observed in ptrans are much larger than the ones 
observed in pi ong , since for the first switching process (B) ptrans jumps from its minimum 
value at 315° to its maximum value at 225°, while pi on g exhibits the same values at 315° 
and 225°. Due to the larger resistivity difference in p trans in the following we will restrict 
our discussion to this quantity. Note that we nevertheless observe a consistent behavior for 

Plong • 

Figure [HJ^c), (a) shows the dependence of the first (B) and the second (D) switching field 
on Vp for the same magnetic field orientation a = 350°. For positive (negative) V p the 
first switching field shifts to lower (higher) and the second switching field to higher (lower) 
absolute values of the magnetic field. This observation can be explained via a slight tilt of 
the relative orientation of the magnetic easy axes. A positive V p induces a uniaxial magnetic 
anisotropy with a magnetic hard axis along j leading to a switching angle < 90° for the first 
and > 90° for the second switching process. The smaller the switching angle, the smaller 
is the energy cost for the domain wall formation around a new magnetic domain nucleus. 
The latter has to be accounted for by the energy gain of the magnetization switching into 
the energetically more favorable minimum in free energy. Therefore, a smaller (larger) 
angle between two magnetic easy axes leads to a shift of the corresponding switching field 
to smaller (larger) absolute values of the magnetic field, in agreement with experimental 
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FIG. 8: (color online) (a),(c) Enlarged sections of the ptrans(W)-£0 curves measured at 5 K for the 
external magnetic field H aligned at a = 350° showing the piezo- voltage dependence of the first 
(c) and second (a) switching field. (b),(d) For a = 260° the order of switching fields depending on 
the piezo voltage is inverted. The arrows indicate the direction of the field sweeps. 

observation. For the external magnetic field aligned along a = 350° — 90° = 260° the order 
of the switching fields is inverted compared to a = 350° [Figure E(b),(d)]. This is due to the 
fact that the order of the switching angles is inverted, too, corroborating our analysis. 



VII. ADMR AND PIEZO-STRAIN AT 50 K 



Due to the strong temperature dependence of magnetic anisotropy in GaMnAs (Sec. [U and 
IXl) the magnetic anisotropy landscape can be adjusted to a regime in which piezo-voltage 



induced strain qualitatively alteres magnetic anisotropy [50|. This becomes obvious from 
the {piong(tt), Ptrans(tt)} ADMR measurements at 50 K (Fig. [9]). For /i H = 100 mT, the 
magnetization orientation is still predominantly determined by the orientation of the external 
magnetic field a. Therefore, at /j,qH = 100 mT the {pion g (a), Ptrans(«)} curves at different 
piezo voltages only exhibit minor changes. However, for fi H = 10 mT, the magnetization 
orientation and thus also the magnetic anisotropy is strongly influenced by the piezo-induced 
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FIG. 9: (color online) Longitudinal resistivity pi on g(a) and transverse resistivity ptrans(o) at 50 K 
for a constant external magnetic field rotated within the film plane at different voltages V v applied 
to the piezoelectric actuator: (a),(b) V p = +200V, (c),(d) V p = 0V, and (e),(f) V p = -200V. 
The black and grey symbols correspond to the ADMR curves obtained at [iqH = 10 mT and 
fioH = 100 mT, respectively. The yellow (poH = 10 mT ) and red (fioH = 100 mT) dashed 
curves are obtained from the iterative fitting procedure of the ADMR curves at different external 
magnetic field strengths using the parameter values given in Fig. [T3j 
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FIG. 10: (color online) Free energy normalized to the saturation magnetization F/M{j3) as a 
function of magnetization orientation (5 for an external magnetic field [IqH = 10 mT applied along 
a = —25° taking into account the magnetic anisotropy parameters derived from the ADMR at 
different magnetic field strengths at 50 K depicted in Fig. [9j The magnetization orientation is 
expected to change by A/3 = 76° when the piezo voltage V p is varied between +200 V and —200 V. 

strain, which can be deduced from the qualitatively different behavior of {pi O ng(o0; Ptrans(«)} 
at different piezo voltages. As shown in Fig. [9](a),(b) for an external magnetic field of 10 mT 
at V p = +200 V, pi on g and ptrans exhibit abrupt changes for an angle of the magnetic field 
a ~ 0° between h and j, indicating a nearby magnetic hard axis. The smooth changes in 
resistivity around a ~ 90° indicate a nearby magnetic easy axis. In contrast, the curvature 
for Vp = —200 V [Fig. 0(e), (f)] evolves approximately vice versa, indicating a rotation of the 
magnetic easy axis by almost 90°. In a more thorough analysis the anisotropy parameters 
are determined independently for each piezo voltage via the iterative fitting procedure of the 
ADMR curves at different external magnetic field strengths as already described above for 
the measurements at 5 K (Sec. |V|). In this way we are able to quantitatively determine the 
piezo-voltage induced change in the free energy surface at T = 50 K. This analysis yields a 
rotation of the magnetic easy axis by about 70°. 
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VIII. REVERSIBLE PIEZO- VOLTAGE CONTROL OF MAGNETIZATION ORI- 
ENTATION 



As discussed in the preceding section we are able to rotate the orientation of the magnetic 
easy axis by about 70° via application of different piezo-voltages at 50 K. However, since 
at zero external magnetic field the magnetization can decay into magnetic domains aligned 
either parallel or antiparallel to the magnetic easy axis, a piezo- voltage control of magnetiza- 
tion orientation is not granted. To prevent this issue of domain formation, we apply a small 
external magnetic field fixed at a = —25°, which lifts the degeneracy between antiparallel 
and parallel alignments of m relative to the magnetic easy axis. In this way only one global 
free energy minimum prevails, which is oriented along f3 = —80° for V p = +200 V and along 
P = -5° for Vp = -200 V for /i # = 10 mT (Fig. HDD- 

We now take advantage of the direct correspondence between {pi on g, Ptrans} and the mag- 
netization orientation to in situ monitor the evolution of magnetization orientation as a 
function of V p . A measurement of piong(Vp) and ptrans(^p) allows to recalculate the corre- 
sponding magnetization orientation P(V P ) via inversion of Eq. (T5]) and ([3]). For details re- 
garding the inversion we refer to Ref. 50j. Since the data obtained for fi H = 10 mT applied 



along a = —25° have already been discussed in detail in |50J, in the following we will focus 
on the changes observed for different external magnetic field strengths. Figure [TT1(a).(b) 
[(c), (d)] shows piong(Vp) and Ptrans(^p) at T = 50 K for a constant field fi H = 10 mT 
[HoH = 100 mT] applied along a = —25°. The data consist of three full voltage cycles 
(V p = +200 V— > —200 V— > +200 V) in 5 V steps, evidencing the good reproducibility 
of the measurement. The hysteresis is due to the hysteretic expansion/contraction of the 
piezoelectric actuator with V p [Fig. IU(a)]. 

Figure [121 shows the magnetization orientation /3(V P ) as a function of piezo voltage V p at 
50 K obtained from the pi on g data (black open squares) and the ptrans data (red full circles) 
for different external magnetic fields strengths fioH applied along a = —25°. The abrupt 
jumps in ptrans at V p ~ —75 V and V p ~ V for fioH = 10 mT and 20 mT are an artefact of 
the arcsine transformation caused by the fact that the minimal value of ptrans(V p ~ V) = 
—0.07 mficm in Fig. [TUb) is somewhat larger than the value \p-j = —0.1 mOcm in Eq. ([3]) 
determined from ADMR traces at 50 K (cf. Fig. [T3|) . The magnetization orientations 
P(V P ) obtained from pi ong and ptrans consequently are consistent and demonstrate that for 
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FIG. 11: (color online) Evolution of (a) pi ng(V p ) and (b) ptrans(Vp) as a function of the variation 
of the piezo voltage V p between +200 V and —200 V for hqH = 10 mT applied along a = —25° at 
50 K. The solid green line in (a) displays the elongation dependence of the resistivity parameter po 
in /?i on g [compare Eq. (|2j)] as derived from the corresponding ADMR curves, which is subtracted 
from /9i OI1 g for the magnetization orientation determination. (c),(d) Corresponding measurements 
for pqH = 100 mT. In all panels the orange open circles at V p = —200 V, V, and +200 V 
display the resistivity values obtained from the ADMR measurements at a = —25° (cf. Fig. [9]), 
demonstrating the good reproducibility. 



fioH = 10 mT we indeed are able to continuously and reversibly adjust the magnetization 
orientation at will within about 70° solely via application of a voltage to the piezoelectric 



actuator 



50| | . With increasing magnetic field strength, the Zeeman contribution in Eq. (pQ) 



increasingly dominates over the magnetoelastic contribution to the free energy. Thus the 
angular range within which the magnetization orientation can be rotated as a function of 
Vp decreases from Act ~ 70° for hqH = 10 mT to Aa ~ 15° for [iqH = 100 mT. 
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FIG. 12: (color online) The magnetization orientation (3 as a function of piezo voltage Vp at 50 K 
for different external magnetic field strengths p^H applied along a = —25° calculated from p\ on g 
via Eq. ([2]) (black open squares) and ptrans via Eq. ([3j) (red full circles). The good agreement 
between the /3 values calculated from pi on g and ptrans corroborates our analysis and demonstrates 
that for hqH = 10 mT the magnetization orientation can be continuously and reversibly rotated 
by about A/? m 70° solely by varying the piezo voltage between +200 V and —200 V, as expected 
from the free energy contours (Fig. [T0|) . For increasing pqH the maximum angle of rotation is 
decreased to A/3 » 15° for p H = 100 mT. 

IX. TEMPERATURE DEPENDENCE OF THE MAGNETIC ANISOTROPY 
AND THE MAGNETOSTRICTIVE CONSTANT 



From the full set of ADMR {pi on g(a), Ptrans(a)} measurements and the corresponding 
analyses for V p G {+200 V, V, —200 V} at different temperatures (Sec. IVIIj) . we obtain the 
temperature dependence of the resistivity parameters po, pi, P3, and pj and the anisotropy 
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FIG. 13: (color online) Temperature dependence of the anisotropy parameters Buo(V p ), -E>oiOi 
and S c || (a), (b) and the resistivity parameters po, pi, p%, and pj (c), (d) as derived from the 
iterative fitting procedure of the ADMR curves at different external magnetic field strengths. At 
temperatures T < 7c/ 2 ~ 40 K the magnetic anisotropy is dominated by the cubic contribution, 
while at higher temperatures T > Tq/2 the piezo- voltage induced uniaxial magnetic anisotropy 
contribution becomes more important. Parameter pq in (c) was derived from measurements at 
V p = +200 V and = 100 mT, the parameters pi, ps, and p-j are found to be independent of 
Vp within error bars. 

parameters -Biio(Vp), -Boio, and B c \\ shown in Fig. [131 Within error bars we did not observe 
a dependence of pi, p%, and pj on V p . The dependence of the resistivity parameters on 
the external magnetic field has recently been investigated by Wu et al. 58] in the range 



0.5 T< p H < 9 T. For the magnetic field range PqH < 1 T of our ADMR measurements 
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FIG. 14: (color online) Temperature dependence of the magnetostrictive constant (full black 



371 ] via a nanoelectromechani- 



squares) compared to the values obtained by Masmanidis et al. 
cal resonator (open blue squares). 



we also did not observe a magnetic field dependence of pi, p3, and p-?, which within error bars 
is in agreement with Ref. 581 ] . The piezo voltage dependence of po at 50 K is exemplarily 
shown in Fig. [TlTa).(c). 

As typically observed for GaMnAs, the magnetic anisotropy at low temperatures (T < 
Tc/2 m 40 K) is dominated by the cubic magnetic anisotropy contribution, while at higher 
temperatures (T > Tc/2) the uniaxial magnetic anisotropy contributions gain impor- 
tance 24]]. The strong increase of the magnetoelastic contribution T?no(±200 V) in the 
temperature range 25 K < T < 50 K is due to the increase of absolute piezo elongation with 
temperature (Fig. H]). 

In the following we use the magnetic anisotropy constants B c \\, B 010 and B no as derived 
from the analysis of the magnetotransport measurements (Fig. [TBI) and the piezo-induced 
lattice distortion (Fig. H]) to derive the magnetostrictive constant Am. In contrast to -B 110 , 
this quantity is a property of the GaMnAs film and not dependent on the particular method 
used to apply stress. As detailed in the Appendix, the magnetostrictive constant Am(T) 
can be calculated from 

2AB 1W (T)M(T) 



Alll(T) 3c44[A eii (T)-Ae tt (T)]' (4) 
using the temperature dependence of the piezo-induced lattice distortions Aeu = eu(V p = 
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+200 V) - e u (V p = -200 V), i £ {j,t} (Fig. H, of AB 110 (T) = B 1W (T,V P = +200 V) - 
Bno(T, V p = —200 V) [Fig. [T3l(a)]. and of the magnetization M{T) determined via super- 
conducting quantum interference device (SQUID) magnetometry at a fixed magnetic field 
strength of ^iqH = 100 mT. For C44, we use the elastic modulus of bulk GaAs. Since it only 
varies slightly with temperature, a fixed value of C44 = 59.9 GPa 59j is used in all calcula- 
tions. The temperature dependence observed is summarized in Fig. [14] and compared to the 



results of Masmanidis et al. 



37j. Around 10 K to 20 K there is good agreement with the 



values obtained via nanoelectromechanical measurements. Due to the small magnetoelastic 
contribution to the magnetic anisotropy at 5 K, we do not attribute a high significance to 
our data point at this low temperature. The deviation at higher temperatures between our 



data and Ref. 



371 ] can be accounted for by the differences in Tq between the 80 nm thick 



Gao.948Mn .o52As sample in Ref. 
sample with a much higher Tq = 85 K. 



371 ] with Tq = 57 K and our 30 nm thick Gao.955Mn .o45As 



X. PIEZO- VOLTAGE INDUCED NONVOLATILE MAGNETIZATION SWITCH- 
ING 

Finally, we show that the piezo voltage induced strain not only allows to control the mag- 
netic anisotropy, but also can be used to induce an irreversible switching of magnetization. 
The solid black curve in Fig. IToYa) shows the transverse resistivity loop ptrans(/^o-^) at 40 K 
and V p = +200 V, with the external magnetic field oriented along a = 335°. As indicated 
by the black arrows, [IqH was swept from —300 mT to +300 mT and back to —300 mT. The 
resistivity jump in the upsweep at ^iqH m +4 mT is due to a switching of magnetization 
orientation between two minima in the free energy surface. In the free energy diagram for 
HoH « —5 mT shown in Fig. [T5T b) this would correspond to a switching from the minimum 
at f3 ~ 90° into the minimum at ~ 250°. The hysteresis of ptrans(/-*o-fO between up- and 
downsweeps is caused by the magnetization residing in one of these two minima. To induce 
a switching process of magnetization orientation back to (3 ~ 90°, a negative magnetic field 
in excess of the coercive field /io-f^c = — 6 mT has to be applied. As we will show in the 
following, the additional degree of freedom offered by the piezoelectric actuator allows to 
induce a nonreversible magnetization switching already at external magnetic fields below 
the coercive field. 
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FIG. 15: (color online) (a) Transverse resistivity loop ptransif^oH) at 40 K with V p = +200 V and 
the external magnetic field oriented along a = 335°. The open blue circles describe the evolution 
of Ptra.nsifJ'oH) during the measurements sequence A to G described in the text, (b) Piezo-voltage 
induced change in free energy normalized on saturation magnetization F/M{(3) as a function of 
magnetization orientation (3 for hqH = —5 mT. (c) Evolution of /Jtrans during the piezo-voltage 
sweep from V p = +200 V (#i) -> V (B 2 ) -> -200 V (B 3 )^ V (JS 4 )-> +200 V (_B 5 ). 

Starting from the same magnetic preparation (hqH = —300 mT, V p = +200 V) the 
upsweep [blue symbols in Fig. [ToT a)] is stopped in a second experiment at fi H = —5 mT 
(position denoted Bi). Keeping the magnetic field constant at this value, the piezo voltage 
is swept from V p = +200 V to V p = —200 V with the corresponding increase in ptrans 
shown in Fig. dS^c) [Bx(V p = +200 V)-> B 2 (V P = V)^ B 3 (V p = -200 V)]. Sweeping 
the piezo voltage back to V p = +200 V, Ptranst^o-fO remains in this high resistance state 
[B 3 (V P = -200 V)-> B±(V P = V)-> B 5 (V V = +200 V)]. This behavior can be understood 
with the help of the corresponding free energy plots at V p g{+200 V, V, -200 V} in 
Fig. HBTb). These have been derived from ADMR measurements as described in Sec. IVl 
The free energy plots indicate that the magnetization orientation is rotated from (3 f=s 90° 
at V p = +200 V (Bi) over (3 w 135° at V p = V (B 2 ) to w +200° at V p = -200 V (B 3 ). 
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Sweeping V p back to +200 V, due to the potential barrier at (3 « 180°, the magnetization 
does not rotate back in the same way, but evolves into the potential minimum at ~ +250° 
(B 3 — > B± — > B 5 ). Therefore, the sweep in piezo voltage from V p = +200 V to V p = —200 V 
and back to V p = +200 V results in an irreversible magnetization orientation change and the 
corresponding resistivity change observed at HoH = —5 mT. We note that these experimental 
results are closely comparable to those reported in manuscript 5JJ- 

Sweeping [iqH from -5 mT to +6 mT (B^ — > C) after the piezo-loop cycling of V p from 
+200 V to -200 V and back to +200 V, ptrans decreases and approximately follows the resis- 
tivity curve obtained in the conventional magnetotransport measurement for the opposite 
sweep direction in this magnetic field range. The fact that ptrans does not exactly coincide 
with the downsweep curve we attribute to first indications of multidomain effects, which 
begin to become important at these small magnetic fields. In particular, the downsweep 
curve depicted as the full line in Fig. [ToT a) was measured after saturating the magnetization 
at +300 mT, so that multi-domain effects should be less important. As we sweep HqH back 
to -5 mT (C — > B 5 ) and up to +2 mT again (B 5 — »■ D) p tr ans remains on this resistivity 
branch. 

At HqH = +2 mT (D) we performed a second irreversible switching process via cycling the 
piezo-voltage V p again from +200 V to —200 V and to +200 V. In analogy to the switching 
process at —5 mT described above, the magnetization is switched from the minimum at 
(3 ~ 250° into the minimum at f3 « 90° (D — > E). Sweeping fi H to +3 mT (F) after 
this second irreversible piezo-induced switching and back to —6 mT (G), ptrans as expected 
follows the corresponding upsweep branch. These experiments clearly demonstrate that the 
piezo voltage not only allows to control the magnetic anisotropy and thus the magnetization 
orientation, but also to induce an irreversible magnetization switching. 



XI. CONCLUSIONS 



In summary, we have investigated the magnetic properties of a piezoelectric actua- 
tor /ferromagnet hybrid structure. Application of a voltage to the piezoelectric actuator 
results in an additional uniaxial magnetoelastic contribution to magnetic anisotropy. For 
the GaMnAs thin film investigated, the axis of elongation (contraction) corresponds to the 
magnetic hard (easy) axis of this uniaxial magnetic anisotropy. We studied the tempera- 
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ture dependence of the magnetic anisotropy including the magnetoelastic contribution using 
anisotropic magnetoresistance techniques. The temperature dependence of the derived mag- 
netostrictive constant Am is in agreement with the earlier results of Masmanidis et al. 37] . 
At T = 5 K the magnetoelastic term constitutes only a minor contribution compared to the 
dominating cubic term, but we showed that the switching fields of p(fioH) loops are shifted 
by the application of a piezo voltage at this temperature. At 50 K — where the magnetoe- 
lastic term dominates magnetic anisotropy — we demonstrated a continuous, fully reversible 
control of magnetization orientation by about 70°, solely via application of voltage to the 
piezoelectric actuator. This piezo-voltage control of magnetization orientation is directly 
transferable to other ferromagnetic/piezoelectric hybrid structures opening the way to new 
innovative multifunctional device concepts, such as all electrically controlled magnetic mem- 
ory elements. As an example, the piezo-voltage induced irreversible magnetization switching 
demonstrated at T = 40 K constitutes the basic principle of a nonvolatile memory element. 

In this paper, we have demonstrated the manipulation of the magnetization in the plane 
of a dilute magnetic semiconductor. However, also a manipulation of the magnetization from 
in-plane to out-of-plane should be possible. For this, the magnetic anisotropy landscape has 
to be changed by reducing the strong uniaxial anisotropy in growth direction present in 
Gai^Mn^As thin films grown on GaAs. As discussed in Sec. II this can be achieved by 
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221 ]. In the presence of 



growing on a lattice matched Gai-yln^As virtual substrate 
a dominantly cubic magnetic anisotropy with only a small remaining uniaxial anisotropy 
in growth direction, the application of a biaxial tensile strain within the film plane would 
thus allow a switching of the magnetization orientation from in-plane to out-of-plane. Such 
controllable biaxial in-plane strain could be realized e.g. via cementing the GaMnAs film 
onto the head side of the piezoelectric actuator. The controlled application of either uniaxial 
or biaxial stress would allow a full three dimensional control of the magnetization orientation. 

A possible device concept based on this approach is exemplarily illustrated in Fig. [TBI 
consisting of a GaMnAs layer which is strained locally by a piezo- or ferroelectric actuator. 
To allow optimal transfer of the stress from the actuator to the DMS layer, the Gai_ x Mn x As 
thin film should be made self-supporting by removing an AlAs sacrificial layer below it via 
wet-chemical etching. Finally, prestraining of the DMS can be achieved by the Gai^In^As 
virtual substrate discussed above or nanopatterning. To realize efficient strain transducers, 
a direct deposition of the actuator material onto the DMS is required. The fabrication of 
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FIG. 16: (color online) Sketch of a spin valve structure based on strain modulation. The mag- 
netization of the central part of the Hall bar is manipulated by the actuator positioned above a 
Gai-^Mn^As layer prestrained by a self-supporting Gai_ y In y As virtual substrate. In this partic- 
ular realization, the Gai-yln^As buffer with a lattice constant greater than that of the GaMnAs 
to layer leads to an out-of-plane orientation of the magnetization, which is flipped in-plane by a 
biaxial compressive strain induced by the actuator. 

standard piezo- or ferroelectrics, such as PZT or BaTi03, is well established, but typically 
requires temperatures in excess of 400°C 60, |6J|. Therefore, these deposition processes 
cannot be straightforwardly transferred to GaMnAs films, as this ferromagnetic semicon- 



ductor is grown by LT-MBE at substrate temperatures around 250°C or below The 
exposure of GaMnAs to temperatures above 250°C results in Mn segregation and MnAs 
formation and thus to a destruction of the ferromagnetic semiconductor. To circumvent 
this issue, more exotic fabrication routes, e.g. the realization of PZT or BaTiOs via hy- 



drothermal methods 
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631 ] . could be viable approaches. Organic ferroelectrics such as 



poly(vinylidene fluoride/trifluoroethylene) or PVDF/TrFE are an interesting alternative, as 
they can be spin-coated and polymerized at temperatures below 200°C and therefore should 
be compatible with GaMnAs 64l. |65|. 
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APPENDIX: MAGNETOELASTIC CONTRIBUTION TO MAGNETIC 
ANISOTROPY ENERGY 



In this appendix we discuss the influence of the piezo-voltage induced strain on the 
magnetic anisotropy of the film. As described in section II III the (001) oriented sample is 
cemented onto the piezoelectric actuator with the main expansion direction parallel to the 
[110]-direction of GaAs. In the {j, t, n} coordinate system with j| | [110], t| | [110], and n| | [001] 
the strain induced by the piezoelectric actuator in the GaMnAs film is given by the tensor 



e {j,t,n} 



e tt 
e nn 



\ 



(5) 



/ 



with the diagonal components e& describing the strains along the coordinate axes. Note 
that all other components vanish since in the {j, t, n} coordinate system no shear strains are 
present. To derive the piezo-induced change in magnetic anisotropy, we transform e{j,t,n} i n t° 
the {x, y, z} coordinate system with coordinate axes along (100) using the transformation 
matrix 



T 



\ 



(6) 



and its inverted matrix 



1J 




V 



(7) 



I ) 

In this way we obtain the strain tensor in the {x, y, z} coordinate system 



/ i 



e{x, y ,z} = ^e{j,t,n}r 1 
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with the strains e n 



-uu 



\i e ii + e tt) and e z 



z nn and the shear strains e xy 



-yx 



Me- 

2 



e tt ). These are now inserted in the general expression for the magnetoelastic 
contribution to the magnetic anisotropy 661 ] 



magcl 



m , 



ml — 



+ B 2 ( 



(9) 



with rrii denoting the direction cosines of magnetization orientation with respect to the 
cubic axes, the magnetoelastic coupling constants Bi, and describing the strains in the 
{x, y, z} coordinate system. The magnetoelastic coupling constants Bi can be substituted 
with the magnetostrictive constants Ai o and A m using the relations B 1 = f Ai o(ci2 — Cn) 



and B 2 = — 3AmC44 with the elastic moduli q 



Since 



mz + mi 



1, and magnetization is oriented in the film plane (m z 



0), the first term in (Q exhibits no anisotropy and therefore in the following is omitted. 
Furthermore, using m x m y = cos(j3 + 45°) cos(/3 — 45°) = cos 2 (/3) — | and again omitting the 
isotropic part we obtain 

3AmC44(e 7 'j 



magel 



/M 



2M 



6 tt) 2 a T3 2 

cos p = n n o m j- 



(10) 



This expression consitutes a uniaxial magnetic anisotropy contribution along [110] with 
the direction cosine of magnetization orientation relative to the j-axis mj and the corre- 



sponding uniaxial anisotropy parameter Buq 
in the total free energy in Eq. (TXT) . 



3AiiiC44(ejj -—ett) 
2M 



and is used as the third term 
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